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ABSTRACT 

Context. Circumstellar disks are believed to be the birthplace of planets and are expected to dissipate on a timescale of a few Myr. 
' J ' The processes responsible for the removal of the dust and gas will strongly modify the radial distribution of the circumstellar matter 
Ph and consequently the spectral energy distribution. In particular, a young planet will open a gap, resulting in an inner disk dominating 

C/^ the near-IR emission and an outer disk emitting mostly in the far-IR. 

Aims. We analyze a full set of data involving new near-infrared data obtained with the 4-telescope combiner (VLTI/Pionier), new 
mid-infrared interferometric VLTI/Midi data, literature photometric and archival data from VLT/NaCo/SAM to constrain the structure 
of the transition disk around T Cha. 

Q Methods. After a preliminary analysis with a simple geometric model, we used the Mcfost radiative transfer code to simultaneously 

^—^ model the SED and the interferometric observables from raytraced images in the L'-, and W-bands. 

Results. We find that the dust responsible for the strong emission in excess in the near-IR must have a narrow temperature distribution 
with a maximum close to the silicate sublimation temperature. This translates into a narrow inner dusty disk (0.07-0. 1 1 AU), with 
' ' a significant height (H/r ~ 0.2) to increase the geometric surface illuminated by the central star We find that the outer disk starts 
at about 12 AU and is partially resolved by the PioNffiR, Sam, and Midi instruments. We discuss the possibility of a self-shadowed 
inner disk, which can extend to distances of several AU. Finally, we show that the Sam closure phases, interpreted as the signature of 
^ a candidate companion, may actually trace the asymmetry generated by forward scattering by dust grains in the upper layers of the 

outer disk. These observations help constrain the inclination and position angle of the disk to about -1-58° and -70°, respectively. 
Conclusions. The circumstellar environment of T Cha appears to be best described by two disks spatially separated by a large gap. 
The presence of matter (dust or gas) inside the gap is, however, difficult to assess with present-day observations. Our model suggests 
the outer disk contaminates the interferometric signature of any potential companion that could be responsible for the gap opening, 
and such a companion still has to be unambiguously detected. We stress the difficulty to observe point sources in bright massive disks, 
and the consequent need to account for disk asymmetries (e.g. anisotropic scattering) in model-dependent search for companions. 
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1 . Introduction tion is driven by several mechanisms, from viscous transport to 

J, dust growth and sedimentation towards the midplane. These 

Circumstellar disks of young stellar objects are the products of processes will contribute to the dissipation of the circumstel- 

in-falling matter during the process of star formation, and are material. In past decades, thanks to mid-infrared (mid-IR) 

believed to be the birthplace of planetary systems. These mas- spectroscopic campaigns (ISO and Spitzer, for instance), several 

sive gas-rich, dusty disks will dissipate after a few Myr ( Hernan- disks experiencing a transition phase have been discovered (e.g., 

|dez et al.|2007D and will then evolve into gas-poor debris disks, [Strom et al.|1989l [Brown et al.|20'07t [Verhoeff et al.lZoTT) . The 

which contain the leftov ers of planetary forination. As explamed j^ain characteristic of these transition disks is that they display 



in the recent review by [Williams & Ciiiil pOlip , disk evolu- ^ lack of emission in the mid-IR, while they show a strong ex- 

cess in the far-IR. This has been interpreted as a sign of dust 

Send offprint requests to: olofsson@mpia.de clearing in the inner regions, close to the star Several processes 

* Based on Pionier observations collected at the VLTI (European ^an be responsible for the dissipation of the innermost material. 
Southern Observatory, Paranal, Chile) with programs 087.C-0702(B), gravitational influence of a stellar companion (e.g., 

087.C-0709(A), 089.C-0537(A), 083.C-0883(C & D), and 083.C- n a t ^ as v onnov v, , , • 

mnc^A p d\ CoKu Tau 4, .Ireland & Kraus 2008), which can ti'uncate the cir- 

029 j(A & d). ' ■ 
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cumbinary disk. On a timescale of about 10 yr, the disk can 



also be dissipated from inside out by photo-evaporation ( Alexan 
der et al.|2 006 ). The resulting density distribution will produce 



Table 1. Observations log and field-of-view (FoV) for the different in- 
struments. 



a spectral energy distribution (SED) similar to the one observed 
for the transition disk, where emission in excess becomes signifi- 
cant only at long wavelengths. A third possible scenario involves 
the presence of a recently formed planet in the disk, which will 
open a gap both in the gas and dust density distributions by push- 
ing away surrounding material (e.g., Papaloizo u & Lin| [i984). 
The mass of the forming planet will influence the size and depth 
of th e gap (jCrida et a l. 200 6), and the mor phology of the outer 
disk ( [Ayliffe et al.|2012||Pinina et al. 2012). 

T Cha is one of the four transitional disks studied by |Browri| 
et al.| ( |2007| l, and it has a characteristic SED which displays a 



significant lack of emission in the mid-IR. |Brown et al. ( 2007| l 
modeled the circumstellar disk by defining two spatially sepa- 
rated regions, leaving a gap between 0.2 and 15 AU. This gap 
was required to account for the lack of mid-IR emission. In 
|01ofsson et al.| ( 2011| l we presented and modeled VLTI/ Amber 
observations, with which we spatially resolved the inner dusty 
disk. We found this disk to be extremely narrow and close to the 
star (0.13-0.17 AU). Given the small field of view of the inter- 
ferometer (for this dataset, ~ 60 milli-arcsec, mas hereafter) and 
the overall low signal-to-noise ratio (S/N hereafter), the Amber 
data did not provide additional constraints on the outer disk. In 
|Cieza e t"ar] (pOl 1 ), we used Herschel/PACS and Spire observa- 
tions at 70, 160, 250, 350, and 500 /vm to study the structure of 
the outer disk (unresolved by the instruments). Because of the 
steep decrease in flux at far-IR wavelengths, we concluded the 
outer disk is fairly narrow with an outer radius between ~20- 
30 AU, depending on the surface density profile (between -0.75 
and -1, respectively). 

Although the exact geometry of the system is revealed lit- 
tle by little as new observations are performed, they confirm 
there is a significant gap in the density distribution of the disk. 
This structure could be explained by the candidate companion 
detected by Huelamo et al.|pOTT] ), which lies at a projected dis- 
tance of 6.7 AU from the star As the companion pushes away 
surrounding material, it opens a large gap in the circumstellar 
disk, removing the dust content in the intermediate regions of 
the disk. However, models of gap opening by a planetary com- 
panion hardly explain the large size of the gap that observations 
suggest in the case of T Cha ( Zhu et al.|2012[). Thanks to the Pio- 
NiER visitor VLTI instrument ^Le Bouquin et al. 201 l| l, we revisit 
the disk structure. 

In this study, we present new H- and A^-band interferometric 
observations of T Cha (spectral type G8, Alcala et al.|1993 1 ob- 
tained with the PioNiER and the Midi instruments (ILeinert et al.' 



|2003), both installed on the VLTI paguenauereT al. 2010). 
Additionally, we include archival VLT/NaCo/SAM observations 
(Sparse Aperture Masking, IT uthill et al.|2010| ) in our analysis. In 
the following, we first describe in section|2|the observations and 
data processing and analyze the new interferometric data with 
a simple geometrical model (Sect.|3]l. In Sect.|4] we summarize 
the radiative transfer code capabilities and discuss how interfer- 
ometric observables are derived from the code outputs. We then 
present and discuss our best fit model to the inner and outer disks 
(Sect.[5]and[6j respectively), as well as the possible presence of 
a self-shadowed disk (Sect|7]i. We conclude in Sect.|8] 

2. Observations and data processing 

This study makes use of VLTI/Pionier observations (pro- 
grams 087.C-0702(B), 087.C-0709(A), and 089.C-0537(A)), 
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UT1-UT2-UT4 HD 107145 


LR 
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VLTI/MiDi data (program 087.C-0883(C, D)), a rchival measure- 
ments obtained with NaCo/SAM, presented in |Huelamo et al.| 
(201 1 , program 084.C-0755), and VLTI/ Amber data already pre- 
sented in Olofs son et al.| ( |20l l] program 083.C- 0295). Table[T] 
summarizes the observations and the instrumental setup for each 
instrument, as well as their field-of-view (FoV hereafter, which 
we assumed to be 2D Gaussians with diff'erent full width at half 
maximum, FWHM). 

2.1. VLTl/PmmER observations 

Between April 2011 and April 2012, TCha has been observed 
several times with the VLTI/Pionier instrument, with different 
configurations of the Auxiliary Telescopes (AT hereafter). The 
interferometer combines the light from four telescopes in the H- 
band, therefore simultaneously providing six squared visibilities 
(y^ hereafter) and three independent (but four in total) closure 
phases per configuration. We aimed at a coverage of the (m, v) 
plane as complete as possible, to probe different spatial scales 
of the circumstellar environment of T Cha. The minimum and 
maximum projected baselines are 6.3 and 123.8 m (resolution of 
55 and 2.8 mas, respectively). Given the faintness of the star 
{H - 7.9 mag), most of the observations were obtained without 
spectral dispersion (i.e stacking all the //-band light into a single 
detector pixel). The observing run of June 2011, however, was 
performed with the small spectral dispersion (3 channels across 
the //-band) and aimed at measuring precise closure phases to 
search for companions. Combined with the fact that the seeing 
conditions were below average, this explains the large uncertain- 
ties of these data points. Because of the large error bars, we 
spectrally binned the dispersed visibilities of this particular run 
for the rest of the analysis. 

Data reduction was performed using the standard pipeline 
"pndrs" described in Le Bouquin et al. ( 201 l| l. Figure[T] shows 
the calibrated data. From left to right, the upper three panels 
show the (m, v) plane, V^, and closure phases, respectively. 
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Fig. 1. Top row: from left to right: coverage of the (m, v) plan, squared visibilities and closure phase for the Pionier observations. Bottom left 
panel: correlated fluxes for the three Midi observations and the Irs spectrum. Filled red diamonds (baseline UT1~UT4) have been discarded from 
the analysis (see text for details). Bottom right panel: NaCo/SAM squared visibilities as a function of spatial frequency. 



2.2. VLTZ/MiDi observations 

The VLTI/MiDi observations were taken in March 2011 using 
the Unit Telescopes (UT), in PRISM mode (spectral resolution 
of R = 30), and aimed at spatially resolving the distribution of 
polycyclic aromatic hydrocarbons (PAHs) in the circumstellar 
disk of T Cha. Data were processed using scripts dedicated to 
observations of faint targets with Midi (see Burtscher et al.|20T2 
for more details) interfaced with the MIA+EWS softwar^ (Ver- 
sion "2012Jan25"). Bottom left panel of Figure[T]shows the cali- 
brated correlated fluxes of T Cha as a function of wavelength for 
the three baselines and the Spitzer/iRS spectrum (aperture of 3.5" 
for Irs versus ~ 300 mas for the Midi mask). Because TCha is 
relatively faint for the instrument (Fy < 1 Jy in the A^-band), the 
Midi photometric observations were noisy and we prefer not to 
use them to compute visibilities from the coiTelated flux. One 
can note that both UT2-UT3 observations (black open squares 
and black filled circles, respectively), separated by 3 days, ap- 
pear to be consistent with each other. As explained in Burtscher 
et al. (2012 1 uncertainties are given by the quadratic sum of the 
photon noise error delivered by the MIA+EWS software and the 
standard deviation of the transfer function over the night of ob- 
servations. Overall, the final uncertainties are dominated by the 



http : //home . s tr w. leidenuniv. nl/~ j affe/e ws/ 



latter one (between 5 and 10%). However, during the night of the 
201 l-Mar-24 (baseline UT1-UT4), only one Midi calibrated ob- 
servation was performed, meaning we cannot estimate the stan- 
dard deviation of the transfer function. Additionally, the fact that 
the correlated flux is higher for the longest baseline (UT1-UT4) 
is surprising, strongly suggesting an underestimation of the sys- 
tematic eiTors. Because we cannot determine proper uncertain- 
ties for this dataset, we cannot assess if this deviation is real or 
not. 



2.3. NaCo/Sm and VLT// Amber literature data 

The Sam observations of T Cha were performed with VLT/NaCo 
on UT4, on the night of the 2010-Mar-17, with the L 27 objec- 
tive and a 7-hole mask. Data processing is described in |Huelarno| 
|et aLl l l^OlT ), and the calibrated visibilities are shown in the bot- 
tom right panel of Figure[T] The large uncertainties on the visi- 
bilities come from the longer exposure time to increase the signal 
in the originally prioritized closure phase signal. The candidate 
companion reported by Huelamo et al. (2011 1 was inferred from 
the L'-band closure phases. These are sensitive to the depar- 
ture from centro-symmetry of the surface brightness distribution, 
and the authors found a minimum and a maximum of -1 .61 and 
1.69°, respectively. The closure phase signature was attributed to 
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the candidate companion by |Huelamo et al.| ( |201 and modeled 
as such. In the present paper the companion is not included. 

As discussed in Olofsson et al.| ( |201 1[ ), T Cha was observed 
with the VLTI/ Amber instrument, in low-resolution mode (LR), 
during two consecutive nights of June 2009 (18 measure- 
ments on two different triplets, see Table T). D etails for data 



processing can be found in |01ofsson et al.| ( 201 1[ ) 



4. Radiative transfer code and interferometric 
quantities 

4. 1 . The McFosT code 

To have a more comprehensive understanding of the circumstel- 
lar environment of TCha, we used the radiative transfer code 



McFOST ( Pinte et al. 2006 2009 ) to simultaneously reproduce the 



2.4. Photometric literature data 

To construct the SED of T Cha, we used the photometric mea- 
surements, the Spitzer/iRS spectrum as in Olofsson et al.| ( [20TT| l, 
and we included the Herschel/PACS and Spire observations pre- 
sented and modeled in |Cieza et al.| ( 20lT| ). The complete SED 
covers the spectral range between 0.36 /im and 3.3 mm. 



3. Preliminary analysis 



Compared to the VLTI/Amber observations presented in Olofs- 
son et al. ( 201 1 \ the uncertainties of the Pionier in i/-band 
are overall much smaller, however, the decrease in with in- 
creasing spatial frequencies is comparable. This drop in was 
interpreted and modeled as the signature of the inner dusty disk 
in [Olofsson et al.| ( |2011| l. Interestingly, at lower spatial frequen- 
cies (B/A <20xir^d-') the display a "plateau" (V^ < 1) 
which suggests a contribution from an extended component, e.g. 
thermal emission from the disk that is later scattered or photo- 
spheric light scattered by the outer disk (see, Pinte et al. [2008, or 
[TatuUi et al.|201 l| for a similar discussion applied to the case of 
HD 100546). 



3. 1 . Ring model 



We first used a geometric model of an uniform ring (jEisner et al.j 
|2004 ) on the VLTI/Pionier data for a preliminary analysis. The 
uniform ring model is described by three free parameters, its di- 
ameter a (in mas), its position angle PA (0° meaning the semi- 
major axis is oriented north-south), and its inclination / (0° is 
face-on). The star-to-ring flux ratio is derived from the 2MASS 
//-band flux and the Kurucz stellar model obtained from SED 
fitting (see Sect.[4]i. We used a Monte Carlo Markov chain pro- 



SED and the interferometric observables. Raytraced images are 
produced for the H-, L'-, and A^-band and used to compute V^, 
closure phases, and correlated fluxes. In Mcfost a disk is defined 
with the following parameters: the inner and outer radii (^n, rout, 
respectively), an exponent a for the surface density distribution 
(S(r) = I.o(r/ro)"') and the mass of the dust Mdust- The verti- 
cal structure follows a Gaussian distribution (exp[-z^/2//(r)^]), 
where the scale height H{r) is defined with a reference value Hq 
at the reference radius tq. The flaring of the disk is described 
via the exponent /3 (H{r) - Ho{r/ro'f). Finally the dust grain 
properties are fully described by the optical constants (used to 
compute mass absorption coefficients for a given dust compo- 
sition), the minimum and maximum grain sizes (^min and ^max, 
respectively), and a power-law exponent p for the grain size dis- 
tribution (dn{s) oc s^ds, p < 0). Several, spatially separated re- 
gions can be defined in Mcfost. The monochromatic raytraced 
images include the thermal and scattered light contributions. 

To define the stellar photosphere within Mcfost, we used the 
same stellar parameters as in Olofsson et al. ( 201 1 Tgff = 5400 K, 
ii=100pc, Av=1.5). We assumed the geometry of the disk to be 
the same as the one originally described in |01ofsson et al.|f201 l| l, 
with an inner dusty disk close to the star that is responsible for 
the near-IR excess and an outer disk located further away from 
the star that contributes to the far-IR. 

The Mcfost model is hardly comparable to an infinitely flat 
disk as assumed in the case of the geometric model, therefore 
the parameters found using this model (/, a and the PA) for the 
inner disk may be different for the Mcfost model. Given that the 
emission is not fully resolved by Pionier (V^ ~ 0.7), / and PA 
cannot be fully constrained by these measurements. However, 
as discussed later on (see Sect. [6]), the outer disk is in the FoV 
of Sam and its inclination and position angle can be further con- 
strained. In this study, we will assume the inner and outer disks 
to be coplanar 



cedure (emcee package, |Foreman-Mackey et a"L]|2012| ) to com- 4.2. Computing visibilities and closure phases 



pute 120000 models. Probability density distributions for each 
of the three free parameters are then computed from the reduced 
(projection of exp[-x^ /2] along the three dimensions of the 
parameter space). The most probable fit to the data is obtained 
for a diameter fl = 1.9+0.1 mas (0.19±0. 01 AU), a position angle 
of 142° + 4, and an inclination / = 53° ± 2 (the uncertainties cor- 
respond to the width cr of a Gaussian fit to the probability density 
distributions). The inferred inclination of 53° is close to our re- 
sult in jOlofsson et al.| ( |201 1) , and close to the most probable 
inclination found for the outer disk (/ ~60°, [Cieza et al. 2011) , 
suggesting that within the uncertainties, the inner and outer disks 
are coplanar. 

3.2. No detected companion in the H-band 

Additionally, we checked for a possible binary signal in the Pio- 
nier data. We used the procedure described in |Absil et al.| ( |201 T\ . 
We do not detect any signal for a companion, down to a flux ratio 
of 1:100 in the //-band, for separations in the range 10-100 mas 
(1-10 AU). 



With Pionier, we observed TCha using the ATs, which have 
a photometric field of view that can be approximated by a 2D 
Gaussian with a FWHM of about 250 mas. Additionally, even 
in its broad-band mode Pionier is not a monochromatic instru- 
ment, but covers a spectral range within the //-band. One has 
to be careful when computing from raytraced images, espe- 
cially when the photometric FoV is large. We can define the 
photometric flux as 



phot 



Je J A 



I(e, A)dedA, 



and the coherent flux as 



^'coh^ r r 1(0, A)e 
Je J A 



-linOBi 



'^dedA, 



(1) 



(2) 



where I{9) is the source intensity as a function of the centro- 
symmetric coordinate on the sky 0, B and A are the projected 
baseline and wavelength of observations, respectively. Then the 
resulting visibility V will be the ratio Fcoh/^^phot- Under the 
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reasonable assumption that the raytraced image is achromatic 
within the //-band, Eq. ([T|l becomes 



-i 



Fpuot = md9. 



(3) 



However, this simplification cannot be done in Eq. Q because 
of the wavelength dependent term in the exponential. For small 
baselines B, the differences in spatial frequencies within the H- 
band will be small and the effect negligible. However, for longer 
baselines, this effect becomes significant, especially when the 
emission is resolved. To account for this wavelength depen- 
dency, we computed the complex visibilities at different wave- 
lengths (21 points between = 1-28 and /Ij^ax - 2.08 //m) 
assuming the raytraced image at Aq - 1.68 //m to be identical 
within this spectral range (I{0, A) can be approximated by 1(0)). 
The complex visibilities were then averaged before computing 
the y2. 

Concerning Sam and Midi observations, we proceeded in a 
similar way to account for chromaticity of the measurements. 
We additionally made sure to adapt the size of the raytraced im- 
ages to the FoV of the instruments (500 and 300 mas, respec- 
tively). To compute Sam and closure phases, we took 21 
wavelength points between 3.2 and 4.4 jum. To compare our 
model to the Midi data, we computed the correlated fluxes at 
three wavelength points: 9, 10 and 12yum (observed data were 
averaged over a +0.5 fim range) in order to see if we could repro- 
duce the wavelength dependency of the observations. Therefore, 
we generated raytraced images at each monochromatic wave- 
length, computed visibilities at the spatial frequency of the ob- 
servations, over a spectral range of +0.5 yum as discussed in the 
previous paragraph, and multiplied the averaged visibilities by 
the total flux in the images. The spectral range of ±0.5 fim was 
chosen to be the same with which we averaged the observed cor- 
related fluxes. The complete analysis of the Midi data, including 
the PAHs feature, will be presented in Meeus et al. (in prepara- 
tion). 



^ 10- 



Photosphere 
{1300 K) 
{1500 K) 




Fig. 2. Blow-up of the SED of TCha in the near-, mid-IR. Stellar 
photosphere is shown in dotted grey. Two Planck functions at different 
temperatures (T = 1300 and 1500 K) are shown in dashed black and 
solid red, respectively. Both Planck functions are scaled to the Irs flux 
at 6 fim. 



this simple exercise, we only considered a single temperature 
for the dust grains, while they most likely have a distribution 
of temperatures. The inclusion of dust at temperatures colder 
than 1300K will increase the emission in excess at about 15 yum, 
which underlines the complexity of reproducing the sharp "V- 
shaped" dip in the mid-IR. Overall this suggests that a vast ma- 
jority of the dust grains has a single, warm (T > 1300 K) tem- 
perature. This can be explained if the inner disk becomes rapidly 
optically thick (as discussed in Grady et al.|2009 in the case of 
HDl 35344 B). 



5. Modeling the inner disk 

In [Olofsson et al.| ( |201 1} we showed that the inner disk is re- 
sponsible for both the decrease of the and the near-IR excess, 
therefore any attempt to model it has to be led carefully through 
a combined approach on the interferometric data and the SED, 
pending the acquisition of a direct image. In the following we 
underline some of the particular features of the inner disk that 
are inferred from the SED. 

5.1. The near-IR excess 

The near-IR excess around T Cha is already significant in the H- 
band (about 25% of the stellar light), suggesting a large amount 
of warm dust is in the immediate vicinity of the star Figure|2] 
shows a blow-up of the SED in the near- to mid-IR (1-30 fim), 
with the stellar photosphere and the Irs spectrum. Additionally, 
to have a first estimate of the characteristic temperature of the 
dust responsible for the near-IR excess, two Planck functions at 
T = 1300 and 1500 K are shown to mimic thermal emission at 
high temperatures. Both Planck functions are added to the stel- 
lar photospheric model and are scaled to match the Irs flux at 
6 fim. One can see that the near-IR excess is slightly under- and 
over-estimated for temperatures of 1300, 1500 K, respectively. 
Most importantly, both Planck functions slightly over-predict the 
Irs flux in the range 8-15 fim, suggesting a decrease with wave- 
length steeper than a Planck function in the mid-IR regime. For 



5.2. Dust composition 

The Irs spectrum of T Cha does not show any emission features 
associated with silicate grains, especially at lOjum, where most 
spectra of circumstellar disks would show the emission band 
arising from the Si-O stretching vibrations in amorphous silicate 
dust grains (see [Henning, 20 10 for a review of cosmic dust op- 
tical properties). There could be two different ways to suppress 
the 10 fim emission feature: (i) the inner disk does not contain 
any silicate dust grains, (ii) silicate grains are typically larger 
than 5-lOjum. Indeed, once silicate grains have grown to size 
larger than several //m, they no longer display the lOjum emis- 
sion feature. To some extent, this effect can also be mimicked 
by high porosity or carbon inclusion in the silicate matrix (see 
|Voshchinnikov & Henning|2008| for instance). 

First, scenario (i) can be ruled out fairly easily as it would 
be counter-intuitive not to have any silicate dust grains in the in- 
ner disk. Silicate grains have been observed in numerous astro- 
nomical environments: the winds of AGB stars, the interstellar 
medium, protoplanetary disks, and planetary crusts in the solar 
system. To match the near-IR excess of T Cha, dust tempera- 
tures do not necessarily need to be higher than the dust sublima- 
tion temperature (T ~ 1500 K, see Fig. [2]). Therefore, there is no 
immediate need to use only refractory dust grains (i.e., resistant 
to higher temperatures). One should note that we do not (and 
cannot) rule out their presence: there is simply no observational 
evidence that requires their inclusion in the model. 



Article number, page 5 of 



13 



Consequently, we examined in detail the second possibil- 
ity in which the inner disk contains silicate grains with sizes 
larger than 0. 1 jum. We used optical constants from |Draine & 



Table 2. Range of parameters explored for the modeling of the inner 
and outer disks. See Sect.l5.3landl6.1lfor details. 



|Lee| ( |1984) l to represent the emission properties of silicate dust 
grains, and did not include any other dust species. We fixed the 
maximum grain size Smax to 1000//m and varied Smin (always as- 
suming p - -3.5 for the grain size distribution, Dohnanyi 1969V 
For a given dust composition, the maximum temperature reached 
is governed by imi,,; the smaller the grains, the warmer. As dis- 
cussed previously, to obtain a satisfying fit to the near-IR excess, 
dust grains need to be warm enough (1300 < T < 1500 K). 
Therefore, the position of the inner radius of the inner disk is, as 
a first approximation, driven by the minimum grain size. With 
imin ^ 5 fim, the lOfim emission feature is always detected in 
the models. For imin = 5 /urn, the dust has to be located close 
to the star (rjn ~ 0.04 AU). This model results in a poor fit 
to the PioNiER V^, especially at long baselines where the mod- 
eled over-predict the observed squared visibilities. Indeed, as 
shown in Sect.[3] the typical diameter probed by the Pionier data 
is of about 2 mas (ri„=0.1 AU at 100 pc). In other words, dust 
grains must be located at about ~0. 1 AU and have a temperature 
between 1300-1500 K. This condition is required to simultane- 
ously model the SED and the V- and can only be achieved if 
small grains (i < 0.1 //m) are included in the model. 

Similar investigations were led by |Meeus et al. ( 2002] ) who 
modeled the SED of HD 100453. The ISO spectrum for this 
source does not show any lOyum emission feature. The authors 
therefore modeled the spectrum of this source using amorphous 
grains with olivine stoichiometry and carbon grains, both species 
having sizes between 4 and 200 yum. The authors also included 
metallic iron grains with sizes between 0.01 and 2.5 fim. The jus- 
tification for including this dust species is that dust temperatures 
higher than the sublimation temperatures of the silicates and 
carbon grains were mandatory to reproduce the near-IR excess. 
This is not the case for T Cha. We conducted similar tests, trying 
to address the following question: how much silicate can we hide 
in the inner disk before the 10 //m feature becomes detectable? In 
addition to the optical constant of astro silicates from Draine &] 
|Lee (1984), we included amorphous carbon ( |Jager et al.|1998| ir 
We choose amorphous carbon, because its mass absorption co- 
efficients are featureless at all sizes. Grain sizes for carbon were 
set to 0.01 < 5 < 1000 yum, and i^^x - 1000 jum for amorphous 
silicates grains (with p - -3.5 for both dust species). We varied 
■Srain for the astro silicates, as well as the relative mass fraction 
between the two dust species and aimed at reasonable fits to both 
the SED and PioNmR V^. We had to change the disk geometry 
in order to test the robustness of each dust setup. We find that a 
relative mass fraction of 70% of astro silicates, with a minimum 
size of 5 fim leads to a good fit to all the observables. With less 
carbon grains, ri„ has to be too close to the star for the to be 
fitted properly. With smaller i^in for the silicate dust grains, a 
10/im emission feature becomes detectable. 

5.3. Disk parameters 



As discussed in Sect. |5.1[ to match the fast decrease of the flux 
between 2 and 15 pm, dust grains must have a narrow range 
of temperatures, centered around 1300-1500 K. Allowing for 
colder grains in the inner disk would increase the flux at about 
15yum. The best fit model for the inner disk is found within a 
grid of models that aim at reproducing both the SED, for wave- 
lengths between 0.3 and 15yum, and the Pionier visibilities, for 
which we used the entire dataset. For the SED, we used 13 pho- 
tometric points (the fluxes at Irs wavelengths were interpolated 



Parameter 



Range 



Inner disk 



W = rout - nn 
//„(ro = 0.1AU) 

Mdust 

/ 

PA 



[0.06, 0.07, 0.08, 0.09, 0.1, 0.11] AU 
[0.03, 0.04, 0.05,0.07, 0.1] AU 
[0.015, 0.0175, 0.02, 0.0225] 
[0.9, 1,2, 3,4,7.5, 10] x lO"" Mq 
[50, 70]°, 1° steps 
[0, 180]°, 10° steps 



Outer disk 



i/„(ro = 25AU) 
/ 

PA 



[8, 9, 10, 11, 12, 13, 14, 15] AU 
[1.5, ...,2.5]AU, 0.1 AU steps 
[1, 5, 6, 7, 8, 9, 10, 20] x 10"^ Mq 
[50, 70]°, 1° steps 
[-180, 180]°, 1° steps 



at 10 and 13 pm). We searched for the best models by minimiz- 
ing the un-weighted sum of reduced for the SED and Pionier 
V^. The uncertainties quoted for each parameters are derived 
from the marginalized probability distributions (oc exp[-;^f^/2]). 
The search is performed on the following parameters: ri„, rout, 
Hq, and M(just- For each model, raytraced images are computed 
for 21 inclinations, and the are fitted for 19 values of the PA. 
Table|2] summarizes the range of parameters explored within the 
grid. As previous analysis ( |Brown et al.||2007{ |01ofsson et al.| 
|201 l| l have shown that the inner disk is narrow, we fixed the sur- 
face density to a - -1 and the flaring exponent /? = 1 as they 
have little to no influence on the resulting SED or V^. Given that 
the Midi data are used to check for consistency, they are not in- 
cluded in the search for the best fit model. Given the much better 
(m, v) coverage and S/N of the Pionier data compared to the Am- 
ber data published in Olofsson et al. (2011), the latter ones are 
not included in the;^'^ minimization process. 

We find the inner disk to be best modeled for rin = 0.07 ± 
0.01 AU and rout = 0.11 ± 0.02 AU, values consistent with the 
results from the geometric model and those from [Olofsson et aT] 
( [2011 1, where we used a slightly different dust composition and 
the data had a lower S/N. We could obtain comparable fits to 
the Pionier with a larger width for the inner disk (rout up to 
0.15-0.2 AU), but the flux in the near- to mid-IR then becomes 
redder and slightly over-predicts the flux at 10-15 pm, while the 
flux at shorter wavelengths is under-predicted. The scale height 
is a crucial parameter. To obtain a majority of dust grains at 
high temperatures, with very few grains at temperatures below 
~ 1000 K, one can increase the scale height of the disk. A larger 
geometric cross-section of the inner disk increases the number of 
grains directly exposed to the stellar radiation, which are there- 
fore efficiently heated up. At a reference radius ro of 0. 1 AU, we 
obtained a good match to the data by setting the scale height at 
Ho - 0.02 ± 0.0025 AU. For such a model, we obtained a max- 
imum temperature of Tdust ~ 1510 K, close to the sublimation 
temperature of sub- /im-sized silicate dust grains. A total dust 
mass of about 2^j x 10"" Mq is sufficient to match the SED in 
the near-IR. The best inclination and position angle, for the inner 
disk solely, are found to be 54+^ °, and 120° ± 30°, respectively. 
The latter two values will be revisited when modeling the outer 
disk (Sect|6]l. The relatively large uncertainty on the position an- 
gle comes from the fact that the emission is not fully resolved by 
Pionier (V^ > 0.7). 
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5.4. Discussion 



6. Modeling the outer disk 



TCha is not the only transition disk showing a lack of silicate 
emission features in the mid-IR (e.g., HD 135344 B, SR21). In 
the case of T Cha, the maximum temperature reached in the best- 
fit model is slightly higher than 1500 K. We hypothesize that the 
inner radius of the inner disk is at the sublimation radius. The 
smaller, and hence warmer, silicate dust grains reach the subli- 
mation temperature and only silicate grains larger than several 
yum may survive to these temperatures. Small, (sub-) yum-sized 
carbon dust grains can survive at higher temperatures, and are 
required in our best-fit model to match the near-IR excess at a 
distance of about 0. 1 AU, a distance constrained by the inter- 
ferometric observations. Overall, given the number of free pa- 
rameters of the model, and that inferring the dust composition 
on a non-detection of the silicate emission feature is not ideal, 
we could hardly investigate the dust properties deeper. We as- 
sumed all grains to be non-porous, hard spheres, while porosity 
and possible inclusions of different materials are known to have 
an impact on the emission features in the mid-IR. 

To reproduce the data, we need the height of the inner wall 
(H/r = 0.2) to be much larger than what hydrostatic equilib- 
rium would prescribe {H/r ~ 0.02). This echoes the results of 
|Thi et alT]p011| . The authors used the thermo-chemical code 
Prodimo to study the morphology of the innermost regions of 
disk around more massive Herbig Ae stars. They investigated 
the influence of decoupled gas and dust temperatures. They 
found the scale height of the innermost regions of the disks to in- 
crease significantly compared to models where the temperatures 
of gas and dust are equal. When decoupled, the gas temperatures 
are higher than the temperature of the dust grains resulting in a 
higher inner rim for the gas. Small dust grains will efficiently 
be dragged along with the gas, which will also lead to a higher 
scale height for these grains. This finding supports the obser- 
vational results of |Acke et al.| ( |2009] l who modeled the SEDs of 
Herbig Ae/Be stars and concluded that a disk in hydrodynami- 
cal equilibrium cannot reproduce the observed near-IR emission. 
We suggest that it is also true in the case of T Cha where the in- 
ner rim is more puffed up. The spatial distribution of the smaller 
grains, which are well coupled to the gas, may increase in the 
vertical direction and hence provide a larger geometric cross- 
section to intercept the stellar radiation. 



As discussed in [ Schisano et al.| (|2009| and references 
therein), T Cha displays a strong variability up to 3 mag in the 
V-band. The authors discussed the possibility of dusty clumps 
in the immediate vicinity of the central star, that could partially 
occult it and be responsible for the observed variability. Such 
fast variability has also been commonly observed around other 
young stars (e.g., AA Tau, [Bouvier et al.|2007]l. Based on CoRot 



photometric observations in NGC2264, 'Alencar et al. ( 2010| 



found that a high fraction of the 83 TTauri stars in their sam- 
ple (30 to 40% of stars having inner dusty disks) exhibit photo- 
metric time variability comparable to the AA Tau phenomenon. 
The authors suggested this variability is the consequence of a 
warped inner disk periodically occulting the photosphere. This 
non-axisymmetric inner disk may trace the magnetospheric in- 
teractions between the disk and the star To account for such 
variability, |Alencar et al.| ( [20T0| ) suggested the scale height at the 
inner rim to be of the order of H/r ~ 0.3 (assuming a random 
distribution of inclinations). Our finding of a puffed-up inner rim 
in the innermost regions of T Cha, combined with time variabil- 
ity may trace a non axisymmetric, warped inner disk interacting 
with the stellar magnetosphere. 



6.1. Disk parameters 

As mentioned in Sect. [3] the Pionier present a "plateau" 
at low spatial frequencies, that can be explained by an over- 
resolved, extended contribution. Assuming a 2D Gaussian pro- 
file for the photometric field of view of the AT (with a FWHM 
of about 250 mas, t/* - 100 pc) material out to 12.5 AU in radius 
can be detected (depending on the inclination /). Therefore the 
outer disk is in the field of view of the AT and may contribute to 
the drop of at low spatial frequencies. The same is true for 
the Midi data where the outer disk acts as a source of incoherent 
flux and contributes to the drop in correlated flux (FoV of about 
300 mas). With the setup for the inner disk as described above, 
we investigated the geometry of the outer disk. Based on all 



available data and the modeling results presented in Cieza et al. 
( |2011| l, we searched for the best model for the outer disk within a 
grid of models. To limit the number of free parameters, we fixed 
some of the parameters already explored in Cieza et al. ( 201 l[ l, 
namely, the flaring index /3 = 1.1, a - -1, and rout - 25 AU, 
leaving the following free parameters: the inner radius, the scale 
height, and the dust mass Mjust- As for the inner disk, several in- 
clinations and position angles for the outer disk are tested against 
the observations. One should note that the inner and outer disk 
are assumed to be coplanar (same / and PA) in the analysis. The 
range of values explored in the grid are reported in Table|2] For 
the dust, we used the optical constants of [Draine & L ee ( 1984), 
with 5n,in=0.1yum and 5n,ax=3mm, assuming p=-3.5. We did 
not include amorphous carbon dust grains in the dust content. 
As opposed to the inner disk, where the lack of emission feature 
at lOjum required to consider an additional dust population, for 
the outer disk, there is no observational evidence to support such 
a choice. 

The best model for the outer disk is found by minimizing 
the un-weigthed reduced for the SED (from optical to mm 
wavelengths), the Pionier V , and the Sam and closure phases 
(see Sect. |6.2[ in which we discuss the best-fit values for / and 
PA). Because the inner disk is slightly different than the one 
used in |Cieza et al.| ( |20lT especially the scale height), we found 
a lower value for the inner radius of the outer disk ( 12 ± 2 AU) to 
match the increase in the SED long-wards to ISfim. The scale 
height Hq at a reference radius of r^ = 25 AU is found to be 
of about 2.1 + 0.2 AU, a value close to the one found in |Cieza| 



et al. 



( |2011| l. We obtained a total dust mass of 8 x lO^M© 
(5 X 10"^ < Mdust < 2 X 10""* M0). With the outer disk being 
in the FoV of Pionier at the AT, it introduces a drop in V- of 
about 5% for the lowest spatial frequency in our dataset. Not 
including the outer disk returns almost equal to unity at these 
low spatial frequencies, which would be inconsistent with our 
observations. 



6.2. The outer disk seen in scattered tight 

When observing with Sam the FoV is of about 500 mas (50 AU 
in diameter at 100 pc). Previous studies of TCha, and results 
from the geometric model suggest that the disk is seen at a rel- 
atively high inclination (/ ~ 55°, see Sect. |5.3| l. Therefore the 
wall of the outer disk {roui - 12 AU) is clearly visible in ray- 
traced images, as well as a large portion of the disk. Because 
closure phases trace departures from centro-symmetry of the sur- 
face brightness, we investigated the closure phase signal gener- 
ated by the disk alone. We computed L'-band closure phases for 
baseline triplets of the Sam observations, for different inclina- 
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Fig. 3. Reduced map and contours, from fitting the Sam closure 
phases, as a function of the inclination ; and position angle PA. 



tions (50° < / < 70°) and position angles (-180° < PA < 180°). 
One should note that for PA = 0°, the Eastern part of the disk 
is inclined towards the observer Figure|3] shows the reduced 

map as a function of both parameters. The minimum re- 
duced Xi for the Sam data is 1.9, found for / = 58+^° and 
PA - -70° + 20°. One should be aware that the Pionier data 
provided a measurement for the PA of the inner disk with a 180° 
degeneracy (PA - 120° or -60°). As a comparison, the reduced 
X^ for a single star (i.e., no closure phase signal) is 3.2. From 
Fig. [3] one should note that for approximately the same inclina- 
tion, a local minimum in the Xr rnap is seen for a position angle 
of about PA = 40° (x^ ~ 2.3). Even the good (m, v) coverage of 
the Sam observations (see Fig.[5]and its description later in text) 
is not sufficient to unambiguously distinguish between these two 
solutions, underlining the need for direct imaging observations. 
Figure|4] shows the observed and modeled c losure phases with 
the binary solution found by Huelamo et al. ( |201 1| reduced Xi 
of 1.5) and our best-fit disk model. Each panel shows consecu- 
tive closure phase measurements (observations were repeated 9 
times) for a given holes triplet. Given the larger FoV and lower 
spatial frequencies of the Sam observations compared to the Pi- 
onier instrument, these observations are best suited to constrain 
the parameters of the outer disk. Consequently, we adopt the 
aforementioned values for / and PA, assuming that the inner and 
outer disk are coplanar, a reasonable assumption given the un- 
certainties we find. 

To have a better understanding of the origin of the closure 
phase signal, we computed raytraced images in which the Mc- 
FOST code separates different contributions; thermal emission, 
stellar scattered light, and thermal scattered light. The thermal 
emission from the outer disk is not responsible for the measured 
closure phase signal. The outer disk emits negligible amounts at 
L'-band, as it is too cold. However, scattering by the outer disk 
surface of L'-band photons originating from the photosphere and 
the hotter inner disk produces a significant closure phase signal. 
These observations probe forward scattering, and therefore con- 
strain the location on the sky of the front side of the disk, that 
we find to be in the North-East comer One should note that 
the phase function of scattering depends, among other parame- 
ters, on the grain size. Until now, we used imin = 0.1 fim for 
the silicate dust grains, which is small compared to the wave- 



length of observations (3.8//m). To test the robustness of our 
predictions, we increased imin to 1 fim in the outer disk, so that 
the grains no longer are in the Rayleigh limit. With Ins > A, 
scattering becomes even more anisotropic (e.g.. Fig. 3 of Mul-] 
ders et aL]|2013) ). With this model, the quality of the fit to the 
closure phases is still in good agreement with the observations 
(Yj - 1.9 for a PA of -69°). This suggests that a good fit to the 
closure phases is achieved as long as light is preferentially scat- 
tered in the forward direction. Within the grid of models for the 
outer disk, we could find slightly better fits to the Sam closure 
phases (x^ ~ 1 .8) when increasing the scale height of the outer 
disk from Hq = 2.1 to 2.3 at ro = 25 AU. However, such models 
overpedict the flux at about 30/vm by a factor < 2. Under the 
assumption the Sam observations trace the outer disk, it could be 
indicative of dust settling in the outer disk: if //m-sized grains 
are present in the uppermost layers of the disk, without larger 
grains, the SED would remain comparable but anistropic scat- 
tering would also happen at higher altitudes. Such models with 
larger scale height additionally increase the over-resolved contri- 
bution of the outer disk in the Pionier FoV, decreasing by a few 
percentage points the at short baselines. However, without 
further constraints on the outer disk, investigating dust settling 
is out of the scope of this paper. 

Bottom three panels of Figure[5] show the phases (not clo- 
sure phases) for the best-fit disk model, the original data, and 
the binary solution found by Huelamo et al.| ( |201 1[ left to right 
panels, respectively). The phases plotted in these three figures 
are obtained via a linear combination of the closure phases. The 
relation is obtained by inversion of the closure phase to phase 
matrix (composed of 1, and -1). The unknown parameters of 
the phases (e.g., piston, tip-tilt) are set at zero. The fit to the 
data is always performed on the closure phases, but the repre- 
sentation of the phases helps to better visualize and compare the 
observations to the models. The best PA found for the disk and 
the binary models are quite different (-70° vs. H-78°, respec- 
tively). Intuitively, one may expect that the semi-minor axis of 
disk (along which the asymmetry is generated) should be aligned 
with the PA of the binary, and hence would expect a ±90° dif- 
ference between the two best-fit position angles (the PA of the 
disk tracing its semi-major axis along the North-South direc- 
tion). In principle, the 7 holes mask of the Sam instrument can 
probe 21 diff'erent phases. But what is actually measured by the 
instrument consists of 35 closure phases (15 independent closure 
phases). Closure phases do not directly trace the phases, and we 
conclude that closure phases can be reproduced at the same level 
with different models. The differences in PA is not an inconsis- 
tency between the two models. 

In the detection paper, the candidate companion was not de- 
tected in the /Ts-band. To check if our model is consistent with 
this non-detection, we computed closure phases in the A's-band 
with the same baseline triplets as the L'-band data and the same 
model (iniin = 0.1 fim, PA = -70°) and we found smaller clo- 
sure phases compared to the L'-band. Overall, we find that only 
a couple of baseline triplets (13 out of 35) show an absolute clo- 
sure signal larger than 0.5° (and only one triplet with an absolute 
closure phase greater than 1°). Depending on the observing con- 
ditions and baseline triplets (i.e., projected baseline, which de- 
pends on the time of the observations), our model is consistent 
with a non-detection. 



6.3. The contribution of PAHs 

The only emission feature detected in the Irs spectrum of T Cha 
is associated with PAHs, at 11.3/vm. In an attempt to obtain 
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Fig. 4. Observed Sam closure pliases in blaclc open circles, modeled closure pliases using a binary model (green line, solution found by^ 

|et al.|2 01 1 ) and modeled closure phases from our best-fit dislc model (red line, PA = -70°). Above each panels, the corresponding holes triplet of 
the 7 holes mask is reported. 
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Fig. 5. Phases representation in the («, v) plane, for the best-fit Mcfost model (left panel), the original Sam data (middle) and the binary model 
from |Huelamo et ^pOl 1[ right). The color of the circles give the sign of the phases, while their sizes depend on their absolute values. 



the most complete description of the entire system, we included 
the contribution of PAHs in Mcfost. We used optical constants 
from | Li & Draine| ( |200ip , and considered the PAHs to have a 
single size of 3.55 x lO'yum. Because the emission of the PAHs 
is non-thermal, they do not necessarily need to be close to the 
star to contribute at 11.3 yum, and they were therefore included 
in the outer disk. This choice is supported by the fact that the 
11. 3 jum emission feature is not observed in the Midi correlated 
fluxes, as opposed to the Irs spectrum (Fig.[T]i. However, the star 
must contribute in the UV regime to excite the electronic states 



of the PAHs. We parametrized the UV flux of the central star 
via its UV ex cess (Luv/^* = .01, assuming a spectral shape 

L,i oc A^-^, see 



Woitke et al.|2010 1. To reproduce the flux level of 



the observed 11.3 yum feature, we found a mass fraction of 0.5% 
of the total outer disk dust mass to be enough. Including the 
PAHs additionally has a small, but still noticeable impact on the 
PioNiER at low spatial frequencies. As a source of non-thermal 
emission, PAHs rendered the outer disk slightly brighter even in 
the //-band. Including them in the model produces slightly lower 
y2 for B/A<20x 10"^ rad"' . 
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However, the best-fit to the Sam closure phases is consider- 
ably degraded (reduced = 7.3 compared to 1.9) when includ- 
ing PAHs. We found the best PA to be 118°. The source of 
asymmetry is this time backward scattering, and the back and 
front sides of the disk are reversed compared to the previous so- 
lution (PA = -70°). This finds an explanation in the very small 
size (s - 3.55 X 10"** /im) of the PAHs grains that are included in 
the best-fit model. Such small dust grains scatter almost isotrop- 
ically and therefore the wall of the outer disk becomes a sig- 
nificant source of asymmetric emission. Therefore, the surface 
brightness distribution of the asymmetry becomes clearly differ- 
ent (wall vs. forward scattering in the upper layers of the disk), 
which explains the difference in;^'^. This suggests that the small 
PAHs grains are not located at the inner edge of the outer disk 
(see VerhoefF et al.|20l"0 for a similar discussion applied to the 
case of HD 95881). Because we do not know the exact size and 
spatial distribution of the PAHs in the outer disk (to be further 
addressed in Meeus et al. in preparation), we opt not to include 
them in our best-fit model. 

6.4. Discussion 

The assumption that the Sam observations do, in fact, trace the 
disk around T Cha and not a candidate companion, has some di- 
rect implications on the location of the PAHs. These grains can- 
not be located at the inner edge of the outer disk, otherwise back- 
ward scattering is too efficient. This implies that PAHs should 
be located further out in the disk, and we hypothesize they have 
been transported outwards, in the outer disk's upper layers by 
radiation pressure or stellar winds. 

The mass of the candidate companion is not well constrained 



wise the Herschel/PACS and Spire observations would be over- 
predicted. Because no stellar companion was detected farther 



and an upper limit of 80Mjup was found by Huelamo et al. 
( |201 I I), suggesting a relatively high mass companion. Simula- 
tions of a gap opened by a planet suggests the outer disk can 
become eccentric depending on the planet mass. With a mass 
above ~ 3Mjup, the companion should rapidly push away mate- 
rial that was at the loci of the Lindblad resonances (e.g., 2:1). 
These resonances are known to efficiently damp any eccentric- 
ity and circularize the companion's orbit. If most of the gas has 
been evacuated from these regions, the damping effect becomes 
negligible. Additionally, other resonances, such as the 3:1 (if 
not cleared), will on the contrary excite the companion's eccen- 
tricity and the outer disk will therefore quickly lose its axisym- 
metry. If the outer disk is warped, shows spiral structures (e.g., 
HD 1 35344 B, Mut o et al.|2012| l, or prominent asymmetries in a 
preferred direction (e.g., HD 142527, Rameau et al. 2012 |Casas- 
|sus et al.|2 0l2), such asymmetries may also be responsible for a 
significant signal in the closure phases. 

However, the phase function of scattering is a challenging 
question. As already mentioned, we used spherical, non-porous 
grains which may not lead to a proper description of the scat- 
tered light contribution. An accurate description (and subsequent 
modeling) of anisotropic scattered light, asymmetries or clumps 
in the outer disk requires at least spatially resolved images of 
the outer disk at different wavelengths (e.g., in the near-IR with 
NaCo in ADI mode, and in the sub-mm with Alma). 

Concerning the size of the outer disk, as in |Cieza et al. 
(2011), assuming a surface density with an exponent a = -1 
we find that the solution of a narrow outer disk (rout = 25 AU) 
provides a good match to the SED, especially at Herschel/PACS 
and Spire wavelengths. We could find comparable solutions with 
greater values for the outer radius (300 AU) but the surface den- 
sity then has to be extremely steep (a = -3) compared to the 
minimum mass solar nebula and other observed disks, other- 



out than 20 AU (very deep imaging of T Cha by Vicente et al. 
2011), this result of a narrow disk cannot be explained by an 
outward truncation of the outer disk. 

[Pinilla et al.| ( |2012[ ) studied the dust distribution in disks 
sculpted by planets. Usually, mm-sized particles are subjected 
to a rapid inward drift. However, when a gap-opening planet is 
included in their simulations, a pressure bump is created at the 
outer edge of the gap. Such a pressure bump will stop the dust 
particles as they start experiencing a positive pressure gradient. 
These mm-sized dust particles will therefore be trapped around 
the pressure maximum. As a consequence, a significant peak in 
the density profile will appear at the edge of gap, where grains 
start to pile up because of the planet-disk interaction. The high 
concentration of mm-sized grains at the edge of the disk will 
then dominate the emission in the far-IR, mm-sized grains being 
efficient emitters at these wavelengths. The outermost regions of 
the disk (beyond the pressure maximum) will then mostly con- 
tain smaller grains (< 100 /vm) that contribute less to the far-IR 
emission. This could mimic a narrow outer disk, even though it 
extends farther out with a population of sub-mm dust grains. 

According tO |Pinilla et al. ( 20T2) l, the dust grains at the inner 
edge of the outer disk will eventually collide at high velocities. 
They will fragment and produce smaller dust grains which will 
follow the gas inside the gap (efficient coupling), and may refill 
the inner disk. 



7. Size of the gap: a self-shadowed disk? 

In the best-fit model presented above, the inner disk is extremely 
narrow, and the gap between the inner and outer disk is there- 
fore very wide. Assuming a single planet lies in the gap, theo- 
ries of planet-disk interactions cannot explain our observational 
result. The size of the gap opened by a single planet, in the 
gas profile, can firstly be approximated as twice the Hill radius 
fH = rp(Mp/3M-i,y^^ . Even with a mass as larg e as 80Mjup, at 



a projected distance of 6.7 AU as suggested by Huelamo et al. 
pOTT), this would lead to a gap size of about 3.5 AU (Zhu et al. 
2012). Within the uncertainties on the position angle of the 
planet relative to the disk and the inclination of the entire sys- 
tem, the inner radius of the outer disk is consistent with this 
value. However, the outer radius of the inner disk does not fit 
in this picture. 

We find the inner disk to be optically thick, with an inte- 
grated optical depth in the midplane of t = 2 in the //-band. 
Material can therefore lie in the shadow of the inner disk. To 
address the question of how much mass can be hidden, we in- 
cluded an additional disk in the Mcfost best model, right behind 
the inner disk, to mimic an larger, self-shadowed disk. For the 
dust composition, we also used the optical constants of 'DraineJ 
|& Lee ( |1984| l, with imin = 0.1 fim and s^ax - 1 mm, assuming 
p = -3.5. We ran several models where we increased the outer 
radius between 1 and 5 AU by steps of 1 AU. For each rout value, 
we increased the mass of dust until the fit becomes inconsistent 
with the data. It appears that the SED and the Midi data are 
best suited to constrain the dust mass in the self-shadowed disk. 
Indeed, the Pionier data are sensitive to the inner radius of the 
inner disk, while the SED will immediately reveal if there is any 
emission in excess in the mid-IR. For a given rout, we stopped 
increasing the dust mass once the modeled correlated fluxes are 
no longer within the 3cr uncertainties of the observed correlated 
fluxes. We repeated this exercise for two different scale height 
reference values (Hq - 0.05, and 0.1 AU at ro - 1 AU). Over- 
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Fig. 7. Best-fit Mcfost model (in solid red) to the SED of TCha 
(photometric observations in black, opened circles), including a zoom 
on the 0.5-35 ycca. spectral range. The dashed line is the photosphere. 



Fig. 9. Observed and modeled visibilities (open circles and 

red squares, respectively) for the best fit-model compared to the 
VLTI/Amber observations presented in |01ofsson et al.|p01 1| >. Top and 
bottom panels display K- and //-band data, respectively. 



all, the dust mass that can be hidden in the shadow of the inner 
disk varies from ~ 5 X 10"" up to ~ 2 x 10"'* Mq for r^ut = 1 
and 5 AU, respectively (2.5 and 100 times the mass of the inner 
disk). With the available dataset, we cannot further constrain the 
extent nor the mass of this self-shadowed disk. 



7. 7. The best-fit model 

Figure|6] shows color-inverted raytraced images of the best-fit 
model in the H-, L'-, and A^-band (left to right, respectively), 
with PA - -70°, / - 58°. The FoV of the different instruments 
are shown as dashed circles on each panels: Pionier and Amber 
in the //-band (250 and 60 mas, respectively), Sam in the L'-band 
(500 mas), and Midi for the A^-band (300 mas). A detailed com- 
parison of the best-fit model discussed above with the various 
data sets is presented on Figs|7] to |9] Fig.|7] shows the best-fit 
model to the SED of TCha from optical to mm wavelengths, 
including a zoom on the 0.5-35 fim spectral region. Upper pan- 
els of FigurejH] show the best-fit model, assuming a PA of -70° 
and an inclination of 58°, to the Pionier (left) and closure 
phases (right). Bottom panels of Figure|8] shows the correspond- 
ing Midi correlated fluxes calculated at 9, 10 and 12 yum (left) and 
the Sam (right). Figure|9] shows the Amber data presented in 
Olofsson et al. ( 2011j ) in K- and //-band (top and bottom pan- 
els, respectively) along with the best fit-model for comparison. 
Even though the best-fit model is in decent agreement with the 
data, the modeled visibilities tend to over-predict the observed 
V^. This may be the consequences of the source variability in 
the innermost regions (e.g.,' Espaillat et al.|2011 1. Table|3]sum- 
marizes the disk parameters and dust composition used in the 
best-fit model. Finally, Figure[TO]shows the temperature map for 
the inner and outer disks. 



0.5 



8. Conclusion 

In this study, we presented a detailed Mcfost model that can ac- 
count for all available interferometric and photometric observa- 
tions of the 7 Myr old transitional disk around T Cha. We find the 
geometry of the circumstellar material to be (at least) twofold. 
First, a narrow inner disk must be located close to the star and 
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Fig. 10. Temperature maps for the inner and outer disks (left and right, 
respectively). The dark contours represent [0.6,0.5,0.4,0.3,0.2,0.1] 
times the maximum density, and the white dashed line represent where 
the integrated radial optical depth equals 1. 



is responsible for the near-IR excess as well as the drop in 
at high spatial frequencies. Based on the significant near-IR ex- 
cess, we concluded its scale height to be important, increasing 
its geometric section, so that a large number of the dust grains 
have high temperatures (~ 1500 K). The lack of solid-state emis- 
sion features in the mid-IR Irs spectrum of T Cha suggests that 
yum-sized silicate grains are at their sublimation temperature and 
are evaporating. Small carbon grains (with a higher sublimation 
temperature) may survive in these regions of the disk and help 
reconciling the interferometric observations and the SED. We 
also discussed the possibility of a more extended (up to several 
AU) self-shadowed inner disk. 

Second, we find that the outer disk (possibly narrow) starts 
at around 12 AU, and is responsible for the far-IR excess. Its 
presence in the Pionier FoV slightly decreases the V' at low spa- 
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Fig. 6. Color-inverted raytraced images of the best-fit model, in the H-, L'-, and A'-band (left to right, respectively). The position angle of the 
disk is -70°. The front side of the disk is in the North-East direction, where forward scattering is observed, while the back side and the wall of the 
outer disk are in the South-West direction. On each panels the dashed circles represent the FoV of the different instruments: Pionier and Amber 
for the //-band (large and small circles, respectively), Sam for the L'-band, and Midi for the W-band. 



Table 3. Parameters for the best-fit Mcfost model. Parameters without uncertainties were fixed in the search for the best-fit model. See text for 
details. 
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tial frequencies via over-resolved scattered light emission. The 
small extent of the outer disk, inferred from SED modeling could 
be the consequence of mm-sized grains accumulating in a pres- 
sure maximum at the edge of the gap, caused by a relatively mas- 
sive companion (> 1 Mjup) opening a gap. The outer disk may, 
however, be more extended in the population of < 100/im-sized 
grains. 

Finally, our best-fit model indicates that the disk itself is ca- 
pable of producing a significant closure phase signal, due to for- 
ward scattering. With present-day observations, constraints on 
the outer disk's morphology, and the scattering phase function, 
we suggest that the Sam closure phase measurements trace the 
outer disk in scattered light and provide a good constraint on its 
position angle and inclination. Several key features of the disk 
around TCha suggest a gap is being opened by a planet, that, 
however, still needs to be unambiguously detected. More gener- 
ally, we stress the overall challenge of observing companions in 
massive bright disks using small-baseline interferometric mea- 
surements and model-dependent technics, as asymmetries of the 
disk (e.g., anisotropic scattering) can strongly contaminate the 
planet's signatures. 
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Fig. 8. Upper left: Observed and modeled Pionier (black circles, red squares, respectively) for the best-fit Mcfost model. Upper right: 
same for Pionier closure phases. Bottom left: observed Midi correlated fluxes (grey circles) and modeled correlated fluxes at 9, 10, and 12 pm (red 
squares). Bottom right: observed and modeled Sam V- (black circles, red squares, respectively). 
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